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We study the effective forces acting between colloidal particles trapped at a fluid interface which 
itself is exposed to a pressure field. To this end we apply what we call the "force approach", which 
relies solely on the condition of mechanical equilibrium and turns to be in a certain sense less 
restrictive than the more frequently used "energy approach", which is based on the minimization 
of a free energy functional. The goals are (i) to elucidate the advantages and disadvantages of the 
force approach as compared to the energy approach, and (ii) to disentangle which features of the 
interfacial deformation and of the capillary-induced forces between the particles follow from the 
gross feature of mechanical equilibrium alone, as opposed to features which depend on details of, 
e.g., the interaction of the interface with the particles or the boundaries of the system. First, we 
derive a general stress-tensor formulation of the forces at the interface. On that basis we work out 
a useful analogy with 2D electrostatics in the particular case of small deformations of the interface 
relative to its flat conflguration. We apply this analogy in order to compute the asymptotic decay of 
the effective force between particles trapped at a fluid interface, extending the validity of previous 
results and revealing the advantages and limitations of the force approach compared to the energy 
approach. It follows the application of the force approach to the case of deformations of a non-fiat 
interface. In this context we first compute the deformation of a spherical droplet due to the electric 
field of a charged particle trapped at its surface and conclude that the interparticle capillary force is 
unfikely to explain certain recent experimental observations within such a configuration. We finally 
discuss the application of our approach to a generally curved interface and show as an illustrative 
example that a nonspherical particle deposited on an interface forming a minimal surface is pulled 
to regions of larger curvature. 
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I. INTRODUCTION 

Experimental evidence has been accumulated that electrically charged, /Ltm-sized colloidal particles trapped at 
fluid interfaces can exhibit long-ranged attraction despite their like charges [ESSfllEIESIll- The mechanisms 
leading to this attraction are not yet fully understood. An attraction mediated by the interface deformation was 
proposed 0], in analogy to the capillary force due to the weight of large floating particles [§, However, for the 
particles sizes used in the aforementioned experiments gravity is irrelevant. Instead, one is led to invoke electrostatic 
forces which act on the interface. This feature has triggered investigations of capillary deformation and capillary- 
induced forces beyond the well studied case of an interface simply under the effects of gravity and surface tension 
P 0, [H M m M: [13, H [l9| . These studies have relied almost exclusively on what we shall call the "energy 
approach", which is based on the minimization of a free energy functional obtained as a parametric function of the 
positions of the particles by integrating out the interfacial degrees of freedom, leading to a "potential of mean force" . 
This functional has to include the contribution by the interface itself, by the particles, and by the boundaries of the 
system. Moreover, due to technical challenges the theoretical implementation of this approach is de facto restricted 
to the regime of small interfacial deformations. 

In the following, as an alternative we investigate the "force approach" which follows by directly applying the 
condition of mechanical equilibrium. Our analysis is based on the pressure field n(r) (generated, e.g., by electrostatic 
forces) acting on the interface between two fluid phases. In general, the condition that an arbitrary piece S of this 
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interface is in mechanical equilibrium reads (see, e.g., Ref. [20I) 

/dAe„n + 7* d^etxe„ = 0, (1) 

Js JdS 

where e„ is the local unit vector normal to the interface, is the unit vector tangent to the boundary dS (oriented 
such that et x e„ points towards the exterior of S), dA is the element of the interfacial area, d£ is the element of 
the arclength along the contour dS, and 7 is the (spatially homogeneous) surface tension of the interface. In Eq. ([T]), 
the first term is the so-called bulk force exerted on the piece S by the pressure 11 and the second one is the line force 
exerted on the contour and generated by the surface tension (also called capillary force). This equation is the starting 
point for the subsequent calculations. 

The force approach allows us to obtain new results, to derive previous ones more easily than within the energy 
approach, and to gain additional insight. This approach was employed in Ref. 21 for the special case that gravity is the 
only relevant force and it was shown to give the same results as the energy approach if the deformations with respect 
to a flat interface are small everywhere. For an arbitrary pressure field acting on the interface, in Ref. [l^ we applied 
the force approach in order to obtain the deformation of an otherwise flat interface far from the particles generating it. 
In the following we further illustrate the force approach: In Sec. |TT]we first express the force exerted by the interface 
in terms of a stress-tensor formulation, extending a recent result [2^] to the most general case n(r) 7^ constant. This 
formulation also allows us to establish a useful analogy between two-dimensional electrostatics and the description 
of small capillary deformations of a flat interface. Since this analogy has been already employed by several authors 
in a more or less explicit manner, here we present a thorough discussion addressing not only the issue of interfacial 
deformations, but also that of boundary conditions and of the capillary forces. In Sec. lIIII we exploit the stress-tensor 
formulation and the electrostatic analogy in order to study the interface-mediated effective force between colloidal 
particles trapped at a fluid interface and to provide a detailed comparison with the corresponding results obtained 
within the energy approach. In Sec. IIVI we compute the deformation of a spherical droplet due to the presence of a 
charged particle at its surface, generalizing a corresponding result obtained in Ref. [l^ for a flat interface and correcting 
certain claims in the literature. Finally we discuss the more general case that the unperturbed interface is curved. 
Sec. |V] provides a summary and an outlook. 



II. STRESS-TENSOR FORMULATION AND ELECTROSTATIC ANALOGY 

The capillary force exerted by the interface (second term in Eq. ([T])) can be rewritten as 

7 ® d^ Gf X e„ = - f d£ {et X e„) • T, (2) 
Jas JdS 

which serves to define the stress tensor T(r) := — 7l(r), where l(r G S) is the 2D identity tensor on the tangent 
plane of S at each point r. In these terms, the condition of mechanical equilibrium (Eq. ([1])) takes the form 

/ dA e„n ^ <f di (et X e„) • T. (3) 

Js JdS 

We recall that ej x e„ is a vector tangent to the surface S but normal to the contour dS and pointing outwards. This 
allows one to reinterpret Eq. ([3]) as the definition of the stress tensor T, the flux of which through a closed boundary 
is the bulk force (first term in Eq. ^) acting on the piece of interface enclosed by that boundary. In dyadic notation 
one has (summation over repeated indices is implied) 

T= -71= -TS^'e^e, ^ T^" ^ - ^ (4) 

where {81,82} is a local basis, at each point tangent to the surface S, and gab = ■ ^b is the induced metric (5"'' 
are the contravariant components of this tensor). In this form, we have the same stress tensor as the one derived in 
Ref. [2^ using methods of differential geometry within the energy approach and for the restricted case n(r) — constant. 
We remark, however, that Eq. ([3]) holds for an arbitrary pressure field n(r). 

An analogy with 2D electrostatics emerges by considering small deformations relative to a flat interface^ corre- 
sponding to the generic experimental set-up, i.e., a situation like the one in Fig. 3] but with, e.g., a charged particle 



^ More precisely, the small quantity is the spatial gradient of the deformation (see Eqs. |[5l l6t). 
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P, 0, H, H, 0, H, nil , a nonspherical particle [1, H^I , or a droplet at a nematic interface [1^ , so that the interface is 
deformed by an electric field, by a nonplanar contact line, or by the elastic stress in the nematic phase, respectively. 
To this end, we identify the flat interface with the XY-plane, so that any point of the deformed interface can be 
expressed as r = r|| + e^u(r||) with r|| := xe^ + yey, and (the subscript || will be used to denote quantities evaluated 
at and operators acting in the reference, flat interface) 



(5a) 



etde = (d£|| • V||)r = d^ii +e^(d£|| • V||u), (5b) 

where d^y = e^dx + eydy, dAy = dxdy, and Vy = ex{d/dx) + ey{d/dy). We denote the projection of any piece 
of interface S onto the Xl"-plane as S\\ , and introduce the unit vector n in the XY-plane which is normal to the 
contour dS\\ and points outwards. With this, we expand Eq. ([3|) in terms of the deformation u(r||); to lowest order 
the component in the direction of is linear in u{rn), 



dA||n=-7* d^||n-V||W. (6a) 
S|| Jasn 

The local version of this equality is the linearized Young-Laplace equation: 

7V|m = -n. (6b) 

To lowest order the components of Eq. ([3]) in the XY-plane are quadratic in the deformation: 

dA|| n V||M= i d£|| n-T||, (6c) 

fill JdS« 



where 

T|| :=7 



(V||w)(V||^i)-ilV,[u|2l|| 



(6d) 



is a stress tensor defined in the XY~p\a,ne. We remark that Eq. ([6c)l also implies Eq. (j6bp upon applying Gauss' 
theorem, demonstrating consistency. 

The form of Eq. (|6b[) allows us to identify u with an electrostatic potential (and — Vyu with an electric field), 11 
with a charge density ("capillary charge" [2^). and 7 with a permittivity^. The boundary conditions usually imposed 
on u at a contour Cy have a close electrostatic analogy, too (see Fig. [T] for the notation): 

(i) The potential is given at Cy. <^ The interface is pinned at the contour at a height z{r^^): 

u(ry) = z(ry), ryeCy. (7) 

(ii) The normal component of the electric field is given at Cy. <^ The contact angle 6'(r||) is specified at each point of 
the contour: 

(n- V||)M(r||) =cota(ry), rySCy, (8) 

where the angle Q!(r||), defined in Fig. [1] must be close to tt/2 for reasons of consistency with the approximation of 
small deformations. 

This means that the equivalence with the corresponding electrostatic analogy is exact concerning the relationship 
between the deformation field u(ry) and its sources (i.e., the pressure field 11 and the boundary conditions). The 
analogy carries over almost exactly, too, to the elastic forces in the XF-plane ("lateral capillary forces") arising 
from interfacial deformation: according to Eqs. ([5cl I6d[) . Ty corresponds to Maxwell's stress tensor, the flux of which 
through a closed boundary gives the electric force acting on the enclosed charge. However, the force related to the 
deformation is actually the interfacial stress, which is minus the flux of this tensor (see Eq. ([2])). Therefore, the 



^ There is also a dual magnetostatlc interpretation in terms of magnetic fields created by currents along the Z-direction; the correspon- 
dences are e2M(r||) ^ A(r|| ), Vy X (e^n) B(r|| ), 7 ^ 1/mo, e2n(r|| ) ^ j(r|| ). 
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(a) (b) 




FIG. 1: Generic view of the three-phase contact line of an interface with a solid boundary, (a) Top view: projection >S|| onto 
the Xy-plane of the interface (dashed area) near the projected contact line (dashed line). The unit vector n is normal to 
the projected contact line and is directed towards the exterior of S\\ . (b) Side view of the same configuration within a vertical 
plane containing n. The thin horizontal line is the section of 5|| (i.e., the fiat, unperturbed interface), while the dash-dotted 
line is the section of the actual interface S, which can be approximated locally by the tangent forming an angle 9 (contact 
angle) with the solid boundary. The latter is the full thick line, which looks locally like a straight line, in general inclined by an 
angle u with respect to the vertical direction (dotted line). The height of the interface at contact entering into the boundary 
condition (O is :— 2:(r|| ). The boundary condition ((Sjl expresses the slope of the interface at contact with a :— + cj. 



electrostatic analogy is valid up to a reversal of the forces and the peculiarity arises that "capillary charges" will 
attract (repel) each other if they have equal (different) sign. The origin of this peculiarity is that Eq. ([3]), which in the 
spirit of electrostatics can be reinterpreted as a definition of T, is actually a relationship between bulk and capillary 
forces as two physically different kinds of forces. The actual connection of T with a force is Eq. 

We note that the electrostatic analogy holds wherever the deformation of the interface is small, even if there are 
other regions of the interface where this is not true. Such "nonlinear patches" can be surrounded by contours where 
the deformation is small, so that the values of the field u and of its derivatives at these contours play the role of a 
boundary condition for the "linear patches" . This means that the "nonlinear patches" are replaced by a distribution 
of virtual "capillary charges" inside the corresponding regions. In particular, there is a simple physical meaning 
associated with the total "capillary charge" and the total "capillary dipole" . The "capillary charge" Q of a piece of 
interface bounded by a contour C is given by Gauss' theorem solely in terms of the value of the deformation at the 
contour (see Eq. (pa|) ): 

Q= ~j f d^ii n- V||M. (9) 

The right hand side of this equation is minus the capillary force exerted on the piece of interface in the Z-direction. 
This implies that in terms of the bulk force Fbuik and by virtue of the condition of mechanical equilibrium one has 

Q = e, • Fbuik. (10) 

This holds even if the deformation in the bulk (i.e., inside the contour C) is not small. In the same manner, it can be 
shown (see Appendix E]) that the total "capillary dipole" P (with respect to the origin of the coordinate system) and 
the torque Mbuik (with respect to the same origin) exerted by the bulk force are related according to 

P = X Mbuik. (11) 

The electrostatic analogy provides a transparent visualization of small interfacial deformations and ensuing forces 
in terms of a 2D electrostatic problem. We note that in Ref . [13 such a kind of analogy is also established, but between 
the capillary interaction and the 3D electrostatic problem (DLVO theory, see, e.g., Ref. |28|) for that very same 3D 
geometrical set-up. As a side remark we note that if gravity (or the disjoining pressure due to a substrate) is relevant, 
it contributes a pressure field which depends explicitly on the deformation field: Ilgiav = —ju/X^ with the capillary 



FIG. 2: Configuration of two particles (top view) a lateral distance d apart, trapped at an asymptotically flat interface. The left 
half S of the full interface is bounded by the particle-interface contact line Co (generically noncircular), the oriented projection 
of which onto the XY-plane is Co,|| , and a contour C, the oriented projection C\\ of which consists of a piece of the F-axis and 
a circle of radius L oo. The unit vector n is normal to the contour and directed towards the exterior of 5. The origin of the 
coordinate system is located at the midpoint O, while the coordinates rjj, ifR parametrize the plane with respect the center of 
the right particle. 

length A [l^. This replaces Eq. (|6b[) by a different equation which is formally equivalent to the field equation of 
the Debye-Hiickel theory for dilute electrolytes (see, e.g., Ref. 128*) with A playing the role of the Debye length. This 
suggests that extending the electrostatic analogy to this case is possible, but this task is beyond the scope of the 
present effort. 



The capillary deformation may give rise to an effective attraction between two identical particles trapped at the 
interface which could explain certain corresponding experimental observations. We compute this force using the 
electrostatic analogy derived above. The study of this issue will also clarify the general relationship between the 
force and the energy approaches as well as the respective advantages and disadvantages. We consider the equilibrium 
configuration of two particles at an asymptotically flat interface and fixed to be a distance d apart (see Fig. [J). The 
total capillary force Ftotai acting on the left half (which includes the piece S and the enclosed particle) is given by 



If the particle separation d is large enough, we can assume that the deformation of the interface at (but not necessarily 
inside) the contour C is small, so that the lateral force F|| (= component of Ftotai hi the XY-p\a,ne = Ftotai — (Ftotai • 
ez)ez) is (see Eq. (pc)) ) 



In the limit L oo, the contribution from the circular part of C|| vanishes and only the knowledge of the deformations 
at the straight midline part of C|| is required, for which the electrostatic analogy will hold provided the interfacial 



III. EFFECTIVE INTERPARTICLE FORCES 



Eq. Il]): 




(12) 




(13) 
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deformations are small there. Thus, one can try to estimate the "electric potential" u at the midline as d — > oo by 
a multipole expansion. The deformation field uji created by the right half plane behaves asymptotically like (see 
Appendix [B| 

URirR) - - — In — + - — > — — + AuR. (14) 

Zttj tr Zttj Zs[rR)^ 

Here, Tr — r|| — is the position of a point relative to the center of the right particle (see Fig. [5]), C is a 

fixed length determined by the distant boundary conditions, which set the zero point (undeformed interface) of the 
"electric potential" ur, and Qg is given in terms of the 2^-pole of "capillary charge" associated with the particle and 
the surrounding interfacial deformation; if this is small everywhere, it is 

Qs = j aAr {trY e^^'^« n(,i„gie) {vr) + "g'^) , (15) 

with the corresponding charge associated to the particle (which is defined by the multipole expansion (see 

Eq. (jB3p ) applied to the interfacial deformation at the contact line). The multipole expansion is based on the implicit 
assumption that the main source of the deformation field is concentrated at or near the particles. For that reason Qs is 
given by the "capillary charge" distribution for d = oo, i.e., for the single-particle configuration. The asymptotically 
subdominant term Aur in Eq. (|14p accounts for the corrections to this assumption, i.e., "polarization" effects by 
the second particle and the fact that, even in the single-particle configuration, the pressure field 11 is expected to 
decay smoothly asymptotically far from the particles rather than dropping exactly to zero beyond some distance: if 
n(r|| — > oo) ~ r^" (actually n = 6 in realistic models in the case of electric stresses [H, [2^, [SO, [Ml , and n = 8 in the 
case of nematic stresses [l^l), there is the bound ly < n — 2 (see the sum in Eq. (fHl) and Appendix IB|) . 

In general the leading term in the multipole expansion is determined by the "capillary monopole" Qq. Thus 
asymptotically for d —>■ oo, F|| will be given by the "electric force" between two monopoles: 

with the reversed sign as explained in the previous Section, so that it describes an attraction. Since Qq is given in 
terms of the net bulk force according to Eq. ([TO]) , the 1/d-dependence only arises if there is an external field acting 
on the system. For example, Eq. corresponds to the "flotation force" (at separations d much smaller than the 
capillary length) for which Fbuik is due to gravity. On the contrary, Qq will vanish if the system is mechanically 
isolated (so that any bulk forces acting on the particle are, according to the action-reaction principle, equal — but 
of oppositte sign — to any bulk forces acting on the interface and hence Fbuik = 0). This was the key point of a 
recent controversy as the force in Eq. ()16|) was advocated to explain certain experimental observations [5, 25, 33] while 
missing that mechanical isolation (as purported in the experiments) rules out this force [ill, IHi [Hj ESj ■ 

If Qo = 0, the capillary force is determined by higher-order terms in the multipole expansion ((T4|). Mechanical 
isolation implies the vanishing of the net bulk force and torque, i.e., "capillary monopole" and "capillary dipole" (see 
Eqs. pm [TT|) ). Thus in general, F|| will take the form of a force between quadrupoles, i.e., it is anisotropic and scales 
as 

- ' (17) 



An experimental realization of this case corresponds to nonspherical inert particles, so that 11 = and the interfacial 
deformation is solely due to an undulated contact line (for a recent corresponding experimental study see Ref. 24). 
Concerning the possibility to relate corresponding experimental observations and theoretical descriptions we point 
out the difficulty that Q2 and higher "capillary poles" are known only in terms of the deformation field u, in contrast 
to Qo and Qi, which are given by the directly measurable and independently accessible quantities "bulk force" and 
"bulk torque" , respectively. 

Another experimentally relevant situation of mechanical isolation corresponds to the case that the "capillary 
charges" Qs of all orders s > vanish. For example, for an electrically charged, spherical particle in mechanical 
isolation one has Qo = and Qi = by mechanical isolation, and by symmetry a rotationally invariant interfacial 
deformation in the single-particle configuration, giving Qs — also for any s > 2. In this case Eq. (fTl)) reduces to the 
correction Aur and the computation of F|| requires a specific model for n(r||) and a detailed calculation. In view of 
our present purposes, we qualitatively derive only a bound on how rapidly F|| decays as function of the separation d. 

The correction AuR{rR — > 00) has a contribution ^ (trY'^^ already in the single-particle configuration because the 
rotationally symmetric "charge density" I\-{rR —> 00) (?'i?)~" does not have a compact support (see Appendix IB|) . 
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This provides a contribution to the capillary force F|| which is equal to the force between the "capillary charges" at 
each side of but near the midline because the presence of the second particle breaks the rotational symmetry. This 
force can be estimated to decay like because the net "capillary charge" in the region farther than a distance 

d from one particle is ^ Jr^^>d ^^11 ^ ^ and the force between charges decays ^ d~^. 

Additionally, Aur has genuine two-particle contributions "induced" by the second particle which can be modeled 
by means of "induced capillary charges" Qs{d) depending on the particle separation. Generically the dominant 
term will be a "capillary dipole"'^ giving rise to a correspondingly dipole-dipole force ~ \Qi{d)\'^ /d'^ , where the 
"induced dipole" Qi{d) must decay at least like the inducing field. This is caused by various reasons. If Qi{d) 
arises by a violation of the boundary conditions at the contact line, one has Qi{d) = 0{A.Ufi{d)) ^ d'^"", and the 
force decays by a factor \/d^ faster than the contribution discussed above concerning the violation of rotational 
symmetry. But there can also be deviations from the linear superposition of the pressure on the interface. This 
occurs, e.g., if the interfacial deformation is due to electric fields emanating from the particles, so that li (x and 
n(r||)-[n(singic)(ri?j + n(singie)(rL)] cx £;(,ingic)(rfl)£;(,ingic)(rL). In this case one has Qi(c;) = C'(i^(singio) (rf)) ^ d-''/^. 
Thus we can conclude that the lateral force must decay as function of separation asymptotically at least as 

i^ll = ©(^-^(^n-a.n+S))^ (18) 



depending on the value of n. In any case this force decays more rapidly than the expression given by Eq. ([16)1 . 

It is instructive to compare with the force obtained within the energy approach. The latter consists of finding 
the parametric dependence on d of the free energy for the two-particle configuration. 



Kion = / d^ 



Vp,rU (19) 



where the integral is the contribution by the interface and Vpart collects the direct contribution by the particles p^[l9|. 
(Within the electrostatic analogy, the effect of V^art would be replaced by appropriate boundary conditions on the 
"potential" u at the interface-particle contact lines.) Vaicn{d) plays the role of a "potential of mean force" for the 
particle-particle interaction, giving rise to a corresponding "mean force" 

i^menW- -KLnCd) (20) 

upon integrating out the capillary degrees of freedom within thermal equilibrium. One should keep in mind that 
this approach captures only the mean-fieldlike contribution to the mean force. The capillary wavelike fluctuations of 
the interface around the mean meniscus profile generates additional, Casimir-like contributions to the force 0, H^, 
which we do not consider in the following. 

One can distinguish two cases. First, if there are "permanent capillary charges" {Qo ^ if the system is not 
mechanically isolated, or 9^ for some s > 2 for mechanical isolation, as discussed above), F|| coincides with 



Ftnen; See, e.g., Refs. [Tol.fl3l for a derivation of Eq. (fTB|) or Refs. |^ [2I, [sl, [13 for obtaining Eq. ([TT]) in the context 
of the energy approach with the simplifying assumption that the interface deformation is small everywhere^. The 
reasoning presented here extends, however, this result also to the case that the deformation around the particles is 
not small, requiring this only near the midline between the particles, i.e., asymptotically for d — > 00. Furthermore, 
the electrostatic analogy shows immediately that the capillary forces -F|| are asymptotically pairwise additive in a 
configuration with more than two particles provided they possess a nonvanishing "permanent capillary pole" . 

The second case corresponds to the absence of "permanent capillary charges" as described above. This has been 
thoroughly investigated in Refs. [T^ [Tsl. [l9l within the energy approach, and has led to i^mGn(rf) ~ d"^""/^, which does 
not agree with any of the possible asymptotic decays indicated in Ec^. (jl8[) . In order to understand this discrepancy, 
we recall that by definition (Eq. ^T^ ) Fn represents the net force acting on the subsystem formed by the particle and 
the piece of interface enclosed by the contour indicated in Fig. O The work done by this force upon an infinitesimal 
virtual displacement 6d is not related in any simple manner to the change SVmcn, which according to the definition 
in Eq. (|19p will involve the work done by local forces during the rearrangement of the "capillary charges" inside the 
subsystem, so that in general Fn ^ Fmcn- 



^ The net "capillary monopole" of the halfplane must vanish exactly due to mechanical isolation and symmetry reasons. Although the 
net torque on the whole system also vanishes, in this case symmetry considerations do not exclude that the net torque on one halfplane 
is opposite to the net torque on the opposite halfplane, so that a "capillary dipole" in one halfplane is possible. 

* Reference performs an exhaustive comparison of the two approaches for the special case that gravity is the only source of deformation 
and the interfacial deformation is small everywhere. 
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In conclusion one has = F^cn if the d-dependence of Vmcn(rf) is dominated asymptotically by a multipole 
expansion, i.e., the whole subsystem can be replaced by a set of point "capillary poles" : the degrees of freedom related 
to the internal structure are irrelevant and only the separation d and the orientation of the "capillary poles" matter. 
This is related to the validity of the "superposition approximation" [oj usually employed in the energy approach, which 
consists of approximating the deformation field u by the sum of the deformation fields induced by each particle in the 
single-particle configuration. In Ref. 19 it is shown that this approximation is valid if the system is not mechanically 
isolated because in that case, asymptotically for d — > oo, the interface-mediated effect of one particle on the other 
amounts to shift it — together with its surrounding interface — vertically as a whole, i.e., without probing or affecting 
the "internal structure" of the subsystem "particle plus surrounding interface" . 

Finally, it is clear that both F|| , defined in Eq. p^ . and F^^n, defined in Eq. (|^D|) . include a contribution from 11 and 
thus differ from the force acting only on the colloidal particle (which would be the integral in Eq. (|13l) extended only 
along the particle-interface contact line) . If one is interested in physical situations in thermal equilibrium, i^mcn does 
represent the effective force between the particles, i.e., once the capillary degrees of freedom have been integrated out. 
In dynamical situations out of equilibrium, completely new considerations have to be made concerning, e.g., whether 
the capillary degrees of freedom can be assumed to have relaxed towards thermal equilibrium in the dynamical time 
scale of interest. But this discussion lies beyond the scope of the present analysis. 



In the following we shall discuss some applications of Eq. ([T]) for particles trapped at an interface which in its unper- 
turbed state is curved. In Subsec. llV Al we shall first consider the interfacial deformation induced by a single charged 
particle on an otherwise spherical droplet. This configuration is particularly relevant for the experiment described 
in Ref. [^. As explained in the previous section, mechanical isolation rules out a monopolarlike (i.e., logarithmic) 
deformation if the unperturbed interface is flat. Since there has been recently a controversy whether this conclusion 
is altered by the curvature of the droplet, we shall present a thorough analysis for such systems. In Subsec. flVBI the 
application of the electrostatic analogy to a generally curved interface is illustrated. 



We consider a charged spherical particle of radius R trapped at the interface of a droplet which resides on a plate 
(Fig.[3|). This configuration models the experiment described in Ref.|5|in the absence of gravity. Our goal is to compute 
the deformation of the droplet far from the particle. Compared with the energy approach, the force approach has two 
advantages: (i) The result is more general because we have to assume only that the deformation is small far from the 
particle; the usual linear approximation is not required to hold also near the particle, (ii) The boundary condition 
"mechanical equilibrium of the particle" is incorporated easily irrespective of the details how the particle is attached 
to the interface. It will turn out that the implementation of this condition has been the source of mistakes in the 
literature. 

We apply Eq. ([T]) to the piece 5(5') of the curved interface bounded on one side by the particle-interface 
contact line Cq and on the other side by a circle given by the constant latitude 5', C(^) :— {ip — "if < ipi}, so 
that dS = Co U C(5'), and we assume that the particle is located at the apex opposite to the plate (Fig. [3]). The 
unperturbed state corresponds to an uncharged particle which does not exert a force on the interface, so that the 
equilibrium shape of the interface is spherical. In the presence of electric charges, the interface will deform. If the 
particle stays at the upper apex, also the deformed interface exhibits axial symmetry. We rewrite Eq. ([1} in three steps. 

(i) The pressure splits into 



Here, i?d is the radius of the unperturbed, spherical droplet, and 2"f /R^ is the pressure jump across the interface in 
the unperturbed state. The dimensionless constant /i accounts for the change in hydrostatic pressure due to enforcing 
the condition of constant droplet volume in the presence of interface deformations, nci(r) is the pressure field created 
by the electric field emanating from the particle which includes electric stresses and an osmotic pressure due to a 
possible discontinuity of the ion concentrations at the interface (see, e.g., Refs.i38l.l3^: this pressure field follows from 
solving the corresponding electrostatic problem. We write 



IV. NONPLANAR REFERENCE INTERFACE 



A. Particle on a spherical droplet 



n(r) = Ap + nei(r), 



Ap = (1 - fi/2) . 



(21) 




(22) 
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FIG. 3: Charged colloidal particle (radius R) at the interface of a droplet residing on a plate. Positions on the interface are 
parametrized by the polar angle ip and the revolution angle ip. Without the colloidal particle and neglecting gravity the droplet 
has spherical shape with radius -Rd and normal vector e^. ry = Rder{ip,'fi) is a point at the unperturbed, spherical interface 
(dashed line) whereas r — +u{'ip,ip)er is a point at the perturbed, nonspherical interface (full line). The perturbed interface 
intersects the plate for tp — ipi. 



fluid phase 1 


r 



^ 










fluid phase 2 











FIG. 4: Side view of a single electrically neutral and spherical particle at a planar fluid interface in equilibrium. The radius 
ro = RsinO of the circular contact line follows from the equilibrium contact angle 0. If the interface was that of a large droplet 
(radius Rd 3> R, see Fig. |3]), this expression would exhibit correction terms of order R/ Rd. In the same manner, the presence 
of charges would deform the interface and introduce corrections which in first order are proportional to the deformation. 
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where we have introduced the dimensionless electric force acting on the whole interface 5incn (which at the apex has 
a hole carved out), 



en := 7; / dA e„nci. (23) 

27r7ro Js^^^^^ 

Here rg = i?sin0, where 9 is the equilibrium contact angle between the fluid phases and the particle; tq is actually 
the radius of the circular contact line of the uncharged particle at a planar interface (see Fig. d]) . 

(ii) Since the particle is in mechanical equilibrium, the contact line force must be balanced by the hydrostatic force 
Fap and the electric force Fei acting on the particle, as well as by any other force Fcxt of external origin (e.g., an 
optical tweezer pushing or pulling the particle): 

7* d^ et X e„ = Fap - 27r7ro(eoi + eoxt), (24) 

JCo 



where we have introduced the dimensionless electric and external forces acting on the particle: 

^ Fci ^ Foxt (25) 



(iii) Finally, the following identity is a consequence of elementary considerations of hydrostatics because Ap in Eq. ^1 
is spatially constant: 



/ d^ e„ Ap + Fap = e,Az («-) Ap, 



/S(*) 

where ^4^(5') is the circular area in the XY plane bounded by the circle C(^). 
Thus Eq. ([T|) can be rewritten as 



(26) 



* d^ et X e„ = 27r7ro(£cxt + £c\ - £n) - e:,A^i^)Ap + dA e„no: 



(27) 



With the notation introduced in Fig. [31 one has (as in Sec. [Ill the subscript y denotes quantities evaluated at and 
operators acting in the undcformed spherical interface, i.e., in tangent planes of the undeformed interface): 



e„ dA 



X ) dipdip = (1 + (l + e,. - 



dip dip J 



Rd J 



Rd J 



dAn 



etde = (d^ii • V|i)r 



1 + 



Rd 



di« 



3^(d£,i ■ V||ii) 



where d^|| = Rd{e^dip + e^sinipdtp), dA|| = R'^ sin ip dijjdip, and 



Vii = 



d 



d 



Rd dip Rd sin ip dtp ' 



(28a) 



(28b) 



(29) 



With this notation, Az{'^) — 7r[i?d + u(\E')]^ sin^ ^f. Equation (|27p can be simplified under the assumption that 
deviations of the actual droplet shape from the spherical one of radius Rd are small in the distant region R/Rd 
^ < ip < ipi, so that the linearized approximation of the deformation is valid and terms quadratic in the quantities 
u, fj,, and Hoi (which vanish in the unperturbed state) can be omitted p^ . We emphasize that this condition does 
not exclude large deviations within the piece iS(^'), in particular near the particle. Due to rotational symmetry, the 
vectorial Eq. ([7f|) is independent of the angle ip and involves only vectors parallel to . One obtains the following 
ordinary differential equation for the function u(p^): 



sin cos ^ 



du 
d^ 



(£cxt + £ci - £n) + T^M-Rd - u sin * 



Rl 



dip sin-0 cosV' nei(V')- 



(30) 
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FIG. 5: The surface <S'^' (5'^') runs along the fluid interface (full line) and the particle (dot) such that it lies in the fluid phase 
exterior (interior) to the droplet. The surface S^^lf^ {^'^Ite) that part of the plate surface which is in contact with the fluid 
phase exterior (interior) to the droplet. The surface Sl encloses the whole system "particle + fluid phases" at a macroscopic 
distance L ^ oo from the droplet. The surfaces are oriented towards the exterior of the corresponding fluid volume which they 
enclose. 



This expression is actually the first integral of the Young-Laplace equation (Eq. (B7) in Ref. [T^ incorporating the 
boundary condition at the contact line. In terms of the functions P{4') ■= cosi/j, Q{ip) := 1 + cosV' lntan(V'/2) and 



S{i;) := f ds sins nei(s) [P(^)g(s) - F(s)Q(^)], (31) 

the general solution is 



m(V') = AP{'^) + ro e, • (sext + £oi - £n)QW + ^A^^d + ^(V-), (32) 

where /i and the integration constant A can be determined by the boundary condition at the plate, i.e., "0 = V'lj 
and the incompressibility condition of the droplet p^ . In the limit i?d oo at fixed r = R^^jj, one recovers the 
results of an unperturbed, planar interface [iGj . We are particularly interested in the deformation given by Eq. (j32p 
in the intermediate range R/Rd ^ '0 ^ 1) from the particle and from the plate. In this range, -P(V') 1, 

Q{tp) ^ Inip + const., whereas 1101(0) will decay in general as implying^ 5(0) ~ 0^^" if n > 2. Thus due to 

the behavior of (3(0) there may be a logarithmically varying asymptotic deformation with the amplitude given by 
ro ez ■ (soxt + £ei - £n)- 

We now consider the situation that there is no external field acting on the system. In this case £ext = and the 
electric forces £ei and en are only due to the charge of the particle. The value of Sei — £n can be obtained by adapting 
the reasoning of Ref. [ij to the droplet geometry. The stress tensor in the fluid phase exterior to the droplet is given 
as — pil + t[,|', where pil is the homogeneous, isotropic stress tensor far from the droplet and T^J'' is Maxwell's 
stress tensor due to the electric field (modified to include the possible osmotic pressure by mobile charges, see, e.g., 

Ref. dH). In the same manner, — (pi + Ap)l + T^^'' is the stress tensor in the interior of the droplet, where Ap is 
given by Eq. (PT]) . With the notations introduced in Fig. [SI the following equations hold: 



L 



.U5i,u5(i) 



dA • T^l^ = 0, (33a) 



^ One finds n = 6 in a realistic model assuming that the whole system has no net charge 1391 . In the presence of a net charge, there is a 
monopolar electric field far from the particle. This interesting case is beyond the scope of the present analysis, the conclusions of which 
only hold if Tl^i{4') decays sufficiently fast. 
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dA . Tif = 0, 



(33b) 



which express that the net force on the exterior fluid phase and on the interior one, respectively, vanishes in equilibrium. 
(The contribution to the integrals from the constant isotropic pressures pi and pi + Ap is zero.) On the other hand, 
by definition one has 

2^7ro(en - £ei) = - / dA • T^J^ - f dA • T^f . (34) 

J5(i) J5(2) 



Combining this with Eqs. ((33|) leads to 



en - Eel 



27r7ro 



dA-T 



(1) 



dA-T 



(2) 



(35) 



where we have taken into account that the contribution of T^j over the surface Sl vanishes in the limit L — » oo 
because the electric field decays to zero far away from the droplet (i.e., we do not consider the possibility that there 
are external electric fields). That is, — (en — £ei) is actually the (dimcnsionless) electric force acting on the plate. The 
calculation of this integral requires to solve the corresponding electrostatic problem. However, on dimensional grounds 
one can obtain the estimate^ |£n ^ £ci| ^ |en|(''o/-Rd)"~^- For an asymptotically planar interface, the logarithmically 
varying deformation due to nonzero values of |en — £ci| leads to a long-ranged effective attraction (see Sec. which 
is the reason that this mechanism has been invoked to be responsible for the apparent attraction reported in Ref. 
If this was the explanation, the measurements in this experiment would imply a value [l3l | \eu — Sd] ~ 10~^. On the 
other hand, the theoretical estimate yields |en — Ecil ~ 10~^|en| with n = 6, so that the experimental results would 
require |en| ^ 10"^. This large value is unlikely for realistic surface charge densities [Tsl, [l3|. 

This result corrects the suggestion made in Ref. [l3 that jen — eci| ~ (ro/i?d)^, inferred from a not applicable force 
balance condition. Indeed, if the deformation is small also at the contact line, the condition "mechanical equilibrium 
of the particle" can be derived from Eq. ([50]) and with an expansion in terms oi tpo = / Rd ^ ^ leads to 



du 

d^ 



RdGz ■ (Scxt + Scl) + Tpa 



1 



-jiRd - u{il:o) 



(36) 



The second term, which is subdominant in the limit —>■ oo, is missing in Eq. (B8) of Ref. 13. We have cross-checked 
this corrected expression by deriving it also within the energy approach employed in Ref. 13, which turns out to be 
algebraically much more cumbersome. 

In conclusion, there persists a logarithmically varying deformation with an amplitude which is very small in the 
limit i?d ~* oo; this is actually a finite-size effect intrinsic to the geometry of the set-up and absent for an unbounded 
flat interface. However, it has the same physical origin as any logarithmically varying deformation of a flat interface, 
namely that the system "particle -I- fluid interface" cannot be mechanically isolated in the configuration of a droplet 



residing on a solid plate. In the absence of the plate one has £n — £oi = due to Eq. ([55]) because >Sp;^jg and S^^^^^ 
are not there and the logarithmic dependence in the range "0^1 disappears. This conclusion corrects a recent claim 
of the opposite in Ref. kOt the relevant errors of this work are pinpointed in Ref. l4ll . in particular the implementation 
of the boundary condition "mechanical equilibrium of the particle" . (Reference |4l| represents incidentally, within the 
energy approach, a further confirmation of our conclusion above.) To facilitate the comparison of our calculations 
with Ref. |40j we make two remarks: (i) The reasoning and results are independent of the precise functional form 
of the electric pressure nci(r); the considerations in Ref. '49 in this respect are thus irrelevant, (ii) As a boundary 
condition for fixing the droplet Ref. 40 employs, instead of a plate at V' = "01 ^-s used here, a fictional pressure field 
ncom('0) Of cosV' constraining the center of mass of the droplet One can easily check that our general solution 
(Eq. ([32])) includes this special case, as the contribution of Hcom('0) in Eq. ((32l) eliminates the singularity of Q{ip) at 
ijj ~ TT and the solution reduces to the corresponding expression in Ref. i40. Thus none of these two issues affects the 
conclusion concerning the logarithmic dependence. 



r(2) 



This estimate is supported by explicit calculations of 11^1 = ■ (T^J'' — t'J*') ■ e„ for realistic models [30l . l39l |. More precisely, the 
electric force exerted on the interface is actually concentrated in a small region of area ~ Tq around the particle, so that one expects 
Tel ~ (7|£nl/^o)-P'(''/''o)i where F is a dimensionless function of order unity at distances r ~ ro from the particle and decaying ~ (ro/r)" 
at distances r ~ ^ ro. On this basis, Eq. I|35|l provides the quoted estimate. 
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B. Particle on a generally curved interface 

If the particle is trapped at a generally curved interface [reference interface), the electrostatic analogy can still 
be exploited provided there is a clear separation of length scales between the typical radius of curvature of the 
reference interface and the size of the region around the particle where the interfacial deformations are appreciable, 
say, roughly a few times the particle size R. Then, at distances from the particle smaller than R^ one can exploit 
the electrostatic analogy in order to study the small deviations of the interface from a reference plane tangent to the 
reference interface at some fixed point near the particle. 

The deviations from the reference plane are given by the displacement field u(r||) = Urci{^\\ ) + Su{r^^), where Mi-cf(r||) 
is the deformation of the reference interface and Su(r\\) is the additional deformation brought about by the presence 
of the particle. Correspondingly, the pressure field can be written as 11 = Hicf + SH. If ^ is a distance from the 
particle beyond which the linearized theory holds (i.e., the deformation near the particle need not be small), then 
in the annulus £ < ry < i?d (with a clear separation ^ -C -Rd so that the following dependence can be observed) the 
solution to the field equation can be written as 



Mr,,) = Bolnf +E^P^^i^^+i[A.c--+A:e-^] 

II s=l I II 

27r7i£<.'<fl, " ^ |r||-rj|| 



(37) 



The constants and Bg are determined by the boundary conditions at r — £ and r — R^ and account for "virtual 
capillary charges" outside the annulus. In particular, we have already seen that 271^3^ is the net bulk force on the 
region r < £ in the direction normal to the reference plane (and in addition to the net force in the reference state 

Mrof(r||)). 

Two new issues arise which we have not considered so far: there is an "external boundary" given by the upper 
bound r|| = i?d and there is a nonvanishing reference deformation Mrof(r||). The first issue is not exclusive for a curved 
reference interface and emerges if the interfacial pressure n(r||) does not decay sufficiently fast with the separation 
from the particle, so that in principle one cannot carry out the limit i?d ~> oo. An example has been studied in full 
detail in the previous subsection, where 5H includes a term —fij/R^ (see Eq. (j2ip ) generated by a nonlocal constraint 
(i.e., constant volume of the droplet). In general, the relevance of the boundary conditions at r|| = i?d introduces a 
nonlocal ingredient preventing general statements based just on the electrostatic analogy localized around the particle. 

The second issue implies an "external electric field", — V||Uicf, giving rise to a new phenomenology. To illustrate 
this point, we consider a reference minimal surface, V|Mrcf = (so that Iliof = 0), containing a mechanically isolated 
nonspherical inert particle (so that SIl = 0) with extension R much smaller than the typical curvature radius i?d of 
the reference interface. As we have seen, Qo = Qi = (i-C, no bulk force normal to the reference plane and no torque 
in that plane), but Q2 ^ 0, so that the particle will experience a force and a torque which are given, to leading order 
in the small ratio R/Rd, by the coupling of this quadrupole with the "external field" Urof- According to Eq. (|B5[) . the 
real-valued quadrupole is characterized by the following 2nd rank tensor: 

Q2 = 92 [{exex ~ GyBy) cos 26 + {e^ey + Bye^) sin 26] , (38) 

where (72 > is the amplitude and 9 is the angle which the principal axes of the quadrupole form with the coordinate 
axes. The origin of coordinates, r|| = 0, is the contact point of the reference tangent plane with the interface and must 
be taken at some point in or close to the particle, e.g., its center of mass. In these circumstances, the electrostatic 
analogy provides the force and the torque, respectively, as 



Q2 ■ V||V||Uref 



(39) 

rii=0 



M = -V|| X 



Q2 ■ V||U 



rof 



rii =0 



(40) 



after reversing the sign with respect to the electrostatic expressions, as discussed above. The most general form of 
the traceless Hessian matrix of the reference interface is given by 



V||V||1trcf(r) 



1 



i?d(r) 



[(e^e^, - eyBy) cos 2(/)(r) + (e^ey + eye^) sin 2(j){r)] 



(41) 
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where i?d(r) > is the absolute value of the radius of curvature and (j){r) is the angle between the principal directions 
and the coordinate axes. Without loss of generality one can choose the orientation of the coordinate axes such that 
(j){r^\ = 0) = 0. By inserting Eqs. ([38l l4T|l into the previous expressions for F and M one finally obtains 

F = g2Cos26iV|| ( ^) + 2^2 sin 261 V||0|^ , (42) 



M = -^^^e.. (43) 
i?d(r,|=0) ^ ^ 

Therefore, the nonspherical particle tends to rotate so that = 0, tt, i.e., in order to align its "capillary quadrupole" 
with the principal directions of curvature of the reference, minimal interface. And when aligned like this, it is pulled 
in the direction of increasing curvature l/i?d. This conclusion complements the result found in Ref. l43l . where the 
energy approach has been applied in order to determine the "potential of mean force" of a spherical^ inert particle in 
a minimal surface. In a multiparticle configuration, this "external" force and this "external" torque compete with the 
capillary interaction between the quadrupoles (see Eq. pT|) ). possibly leading to interesting phenomena concerning 
the 2D patterns formed by the particles. 



V. SUMMARY AND OUTLOOK 



We have investigated the force approach for describing colloidal particles trapped at a fluid interface. This approach 
has allowed us to derive a stress-tensor formulation of the interface-mediated elastic forces for an arbitrary pressure 
field n(r) acting on the interface. In this manner we have been able to generalize some of the results of Ref. |23 obtained 
only for a spatially constant pressure field. It is an interesting, open question whether this result is extendible to, 
e.g., membranes, for which bending rigidity as well as surface tension are relevant, and to other, more general cases 
considered in Ref. [13 and involving constitutive parameters beyond surface tension and bending rigidity. Based on 
the stress-tensor formulation we have worked out a detailed analogy between small interfacial deformations and 2D 
electrostatics, encompassing not only the field equation of the deformation but also the elastic forces transmitted by 
the interface. 

We have exploited the electrostatic analogy in order to compute the dominant contribution to the interface-mediated 
force between two particles if they are far apart. This analogy enabled us to clarify the relationship between the energy 
and the force approach and to reveal the advantages and limitations of each. The definition of the effective force F|| , 
which we have employed in the force approach, differs from the effective force F,ncn introduced in the energy approach 
via a "potential of mean force" . However, the difference is asymptotically negligible if the interfacial deformation in 
the single-particle configuration is expressible as a multipole expansion (i.e., via nonvanishing "capillary poles"). For 
example, if the system is not mechanically isolated or the particles are nonspherical, the force approach allows one 
to extend with relative ease the energy-approach result to cases in which the interfacial deformations are not small 
everywhere. Moreover, it justifies asymptotic pairwise additivity of the force in a multiparticle configuration. One 
must bear in mind that none of the two definitions and i^men is the actual force acting only on the colloid, because 
both take into account the force acting on the surrounding interface. This matters for discussing the dynamics of 
trapped particles. However, in thermal equilibrium -Fmcn is the effective force according to which the equilibrium 
state of the particles is determined. In this situation, the energy approach, which provides -Fmcn, is in principle 
advantageous, while the force approach is more powerful (in the sense that it may facilitate or extend the range of 
validity of the calculations) whenever it can be shown that F^^ « -Fmcn- 

For the experimentally interesting system of a particle at the interface of a spherical droplet in contact with a plate 
(Fig. [3]), we find that the presence of the plate breaks mechanical isolation and leads to a logarithmically varying 
interfacial deformation at distances r from the particle in the intermediate range "particle radius <C r <C droplet 
radius" , with the amplitude of the logarithm vanishing as the droplet radius tends to infinity. Our approach has put 
this finite-size effect on a sound basis. Nevertheless, our numerical estimates show that this logarithmically varying 
deformation is very likely too weak in order to explain the apparently long-ranged attraction observed experimentally 
in Ref. for such a system. However, there are still open questions which we have not addressed but which are 



''' Therefore, the "capillary quadrupole" is not permanent, as in the present illustrative example, but rather induced by the "external 
potential" Mrgf(r||) and, according to Ref. |43|. is proportional to the curvature l/i?d- 
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conceivably relevant for this experiment. We have assumed an electrically neutral system; but if there is a net charge, 
e.g., if the colloidal particle is charged but the droplet is not grounded, additional, long-ranged electric fields arise. 
Another interesting question is the loss of rotational symmetry which occurs if the particle is not fixed at the apex 
of the droplet: this might give rise to an additional force (electrostatic or capillary) pushing the particle towards the 
apex which, in a multiparticle configuration, could be misinterpreted as an effective attraction like the one apparently 
observed also for a planar interface. 

Finally, we have discussed briefly the application of the electrostatic analogy and the associated phenomena arising 
when the unperturbed interface is curved. As an illustrative example, we showed that a nonspherical inert particle 
trapped at a minimal surface is pulled to regions in the interface with larger curvature. 



APPENDIX A: TORQUE BALANCE 



If a piece of interface S is in equilibrium, the total torque on this piece must vanish. In this case, the following 
condition must hold (using the same notation as in Eq. H])): 

/ dA(rxe„)n + 7® d£ r x (e* x e„) = 0. (Al) 
Js Jos 

If the deviations from a flat interface are small, one can simplify Eq. (jAip as in Sec. [Til for the force-balance equation. 
To lowest order in the deformation one obtains 

\ dA\\r\\n + j(h d^ii [r||(n ■ V||-u) - un] I X = 0. (A2) 

This equation implies that the expression in curly brackets vanishes because it is a vector orthogonal to (r|| and n 
lie in the XY plane). In the electrostatic analogy the integral over 11 corresponds to the "capillary dipole" P of the 
piece S. This allows one to rewrite Eq. (jA2p as 



V=--fl d£|| [r||(n- V||u)-un], (A3) 

which generalizes Gauss' theorem (Eq. ©) by expressing the dipole of a region only in terms of the values of the 
deformation field and its derivatives at the boundary. On the other hand, via the general equilibrium condition in 
Eq. (jAip . the right hand side of this equality is related to the torque Mbuik due to the bulk force: 

7 * d^ij [r|| (n ■ V||u) - un] = - x Mbuik- (A4) 

The validity of this expression only requires that the deformation is small at the contour dS^ where the linearization 
is performed, but not inside. This proofs Eq. (|lip . 



APPENDIX B: MULTIPOLE EXPANSION IN 2D 



Here we recall briefly some results concerning the multipole expansion in two dimensions. The "potential" u(r) 
created by a "charge" distribution n(r) is given by (C is an arbitrary constant) 

^^(r) = --^ / dA'n(r')ln^--^ = -77^i?e / dA' n(r') In (Bl) 
27r7 J Q 27r7 J Q 

the second equation introduces the complex variable z — r exp [iip) in order to ease the calculations with Re denoting 
the real part. We first consider the case that 11 has a compact support: n(r) = if ?' > i?. The Taylor expansion 

In(z-z') -Inz-^- I -) (B2) 

s— 1 ^ ^ 

is valid in the complex domain \z'\ < \z\. Inserting this expansion into the general expression (|B1|) one obtains 
straightforwardly 
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valid for r > R, with the complex- valued multipolar charges Qs given by Eq. ((T5)) . As can be easily deduced from 
this latter expression, they can be written as Qs = Qs eiip (is9s), where the amplitude is a positive real number 
and 9s G [0, 2tt) is the angle by which the configuration with the charge Qs is to be rotated in order to achieve a 
configuration with the charge Qs- By using the identity rexp(— — r • (e^ — the expansion (jB3|) can be 

rewritten in a more familiar form involving only real-valued quantities: 



i{r) — m — I y 

27r7 r 2Trj ^-^ 



er (?) e. 



sr" 



'Qs, 



(B4) 



where • indicates s scalar products, e,. = r/r, and 



Qs = Re[Qs{ex - iey) rA {e.^ - ie^)] = gs-Re[e"* "(e^ - iSy) yJ. (e^ - iSy)] 

are the real- valued multipolar charges. 

Assume now that 11 ~ r~" as r ^ oo, so that Qs is ill-defined for s > n — 2. Nevertheless, one can still write 



(B5) 



t r = — In — h - — > ■ 

2ttj r 2ttj ^ 



+ Q*s 



-I- AM(r), 



(B6) 



where p is the largest integer such that i' < n — 2. This expression serves to define Au. By using the Taylor 
expansion (|B2[) again, one can write 



/oo 
dA'n(r') ^ 



1 f z 



'Re 



dA'n(r') 



z\<\z'\ 



In 



1 f z 



(B7) 



In this form one can easily check that Am ^ " for r = |z| ^ oo, and the finite multipole expansion (fTil) holds with 
an extra term which is asymptotically indeed subdominant. 
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